Observation of Spatial Fluctuations of the Local Density of States in Random Media 
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We experimentally study spatial fluctuations of the local density of states (LDOS) inside three- 
dimensional random photonic media. The LDOS is probed at many positions inside random photonic 
media by measuring emission rates of a large number of individual fluorescent nanospheres. The 
emission rates are observed to fluctuate spatially, and the variance of the fluctuations increases with 
the scattering strength. The measured variance of the emission rates agrees well with a model that 
takes into account the effect of the nearest scatterer only. 

PACS numbers: 42.25.Dd, 42.25.Fx, 42.S0.Gt, 32.50.+d 



It is well-known that the spontaneous emission rate of 
an excited quantum emitter is not only a property of the 
emitter itself, but also depends on its surroundings on the 
nanoscale 1, 2\. Control of the emission rate of quantum 
emitters has been demonstrated with a wide range of 
systems, such as reflecting interfaces [H, microcavities 3, 
photonic crystals [B|-0], an d plasmonic nanoantennae 8]. 
The effect of the surrounding of the emitter is described 
by the local density of states (LDOS) that counts the 
number of optical modes available for emission at the 
position of the emitter 0, @|. 

In random photonic media that are promising systems 
to observe intriguing Anderson localization of light 
it is an open question how spontaneous emission rates 
are affected by the surroundings. It has theoretically 
been pre dicted that the LDOS exhibits spatial fluctua- 
tions 13-13}. These LDOS fluctuations are determined 
by light scattering near the emitter and the variance of 
these fluctuations is of the order of the scattering strength 
(X/ki), where k is the wave number of light in the medium 
and I is the transport mean free path. In an infinite 
random medium the LDOS fluctuations are essent ially 



equivalent to the Co intensity correlation function 11 1. 
Co is an infinite-range correlation function for waves in 
random media [14|. The variance of LDOS fluctuations 



was also calculated using a single-scattering analytical 
model and was found to be highly sensitive to the level 
of scattering and absorption at the local scale [l2[ , as well 
as to the dipole orientation [13(. No experimental obser- 
vation of the LDOS fluctuations inside random photonic 
media has been reported so far. 

In this Letter, we present an experimental study of 
fluctuations of the LDOS inside random photonic media. 
We use single nanoscale emitters as internal probes of the 
LDOS. Time-resolved fluorescence is recorded for many 
emitters inside random media with scattering strengths 
1/kt up to 0.12. The distribution of the spontaneous 
emission lifetimes is used to test theoretical predictions 



from several different models. 
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FIG. 1: (color) (a) Time-correlated single-photon counting 
setup with modelocked supercontinuum fibre laser, CCD cam- 
era, and single-photon avalanche photodiode (SPAD). (b) 
Scanning electron micrographs (SEM) of random photonic 
media made of Polystyrene (PS) (top) and Zinc Oxide (ZnO) 
(bottom), (c) Fluorescence image of a single fluorescent 
nanosphere with a diameter of 22 nm embedded in a 7.4(5) 
/^m-thick slab of PS random photonic medium at a depth d 
— 3.3 ^m. 

The experimental setup shown in Fig. [Ha) is designed 
to optically probe a single emitter located deep inside a 
random medium, and to measure its time-resolved emis- 
sion. The emitters are optically excited at a wavelength A 
= 543 nm using a modelocked supercontinuum fiber laser 
with a pulse duration of ~ 5 ps and a repetition rate of 20 
MHz. Excitation was performed through a microscope 
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objective with numerical aperture (NA) of 0.63. The 
emission was collected near A = 620 nm using a dichroic 
beamsplitter (Semrock FF593-Di02-25x36), a cutoff fil- 
ter (Schott RG610), and a 14-nm band pass filter (Sem- 
rock FF01-620/14-25). The location of a single emitter 
is identified using a CCD camera and its time-resolved 
emission is measured by a single-photon avalanche pho- 
todiode with a time resolution of 50 ps. 

We studied two types of random photonic media: slabs 
made from polystyrene (PS) and zinc oxide (ZnO) with 
thicknesses between L = 4.0 and 16.5 /im, see Fig. 
[ljb). To probe the emission rates and thus the LDOS, 
we doped the samples with fluorescent polystyrene 
nanospheres with a diameter of 22(2) nm, an emission 
peak at 612 nm, and a quantum efficiency of 96(4) % 
(Duke Scientific red fluorescent nanospheres). The flu- 
orescent nanospheres are much smaller than the wave- 
length of light and contain about thirty dye molecules 
which all sense the LDOS at essentially the same position. 
Since the orientation of the molecules is random and the 
LDOS is orientation dependent 13, [l5j], dye molecules 



inside one fluorescent nanosphere emit at different rates. 
The fluorescent nanospheres show no photoblinking and 
are unaffected by the chemical environment since the dye 
molecules are protected from ambient oxygen and rigidly 
held in a polymer matrix. To fabricate PS random pho- 
tonic media, a polydisperse suspension of spheres (Duke 
Scientific) doped with fluorescent nanospheres was de- 
posited on a standard microscope cover slide, spread 
uniformly and allowed to dry. This inward-growing self- 
assembly method was previously demonstrated for pho- 
tonic crystal fabrication [l6| . By preparing a polydis- 
perse mixture of different spheres, crystallization was pre- 
vented. To fabricate ZnO random photonic media, a sus- 
pension of pigment (Aldrich ZnO < 1 /im, 99 %) doped 
with fluorescent nanospheres was sprayed on a cover slide 
[l7| . The random photonic media have a low areal den- 
sity of fluorescent nanospheres ranging from 0.03 to 0.12 
/im -2 . As non-scattering reference samples (with 1/kl 
— 0), we prepared fluorescent nanospheres in transpar- 
ent polyvinyl alcohol polymer layers that were spincoated 
on cover slides and covered by Poly (methyl methacry- 
late) index matching layers jl8j | . The transport mean 
free paths for PS and ZnO random photonic media were 
determined by total transmission to be I = 1.7(3) /nm 
and I = 0.82(12) /im at A = 620 nm, respectively, yield- 
ing scatte ring strengths (1/kl) of 0.06 and 0.12, respec- 
tively [17Hl9j ]. For time- resolved fluorescence, we mea- 
sured a total of 96 single fluorescent nanospheres, namely 
12 nanospheres in each of four PS and four ZnO sam- 
ples. We also collected data from 12 single fluorescent 
nanospheres in the reference polymer layer. 

A fluorescence image of a single fluorescent nanosphere 
embedded inside a PS random photonic medium is shown 
in Fig. HJc). The image shows a diffuse spot of the flu- 
orescent light. The background fluorescence of the host 



is negligible compared to the emission of the nanosphere. 
We determined the depth of the nanosphere by modeling 
the diffuse spot of fluorescent light with diffusion theory 
17 , lioj ]. We conservatively estimate a depth accuracy 
of one mean free path. In the particular case of Fig. [2] 
we obtain d = 3.3 ± 1.7/im. For our measurements, we 
selected only single fluorescent nanospheres that were iso- 
lated from other fluorescent nanospheres and situated at 
a depth of approximately half the sample thickness [3] ■ 
In Fig. [H[a), we show the time-resolved fluorescence 
from two fluorescent nanospheres embedded in a PS ran- 
dom photonic medium at the same depth. The two 
nanospheres at different positions yield two different de- 
cay curves demonstrating the spatial fluctuation of the 
LDOS. In Fig.[2][b), the decay curve of a single fluorescent 
nanosphere in a non-scattering polymer layer is seen to be 
purely exponential. The time-dependent fluorescence of 
12 single fluorescent nanospheres inside a polymer layer 
all show exponential decay with the same time constant 
of 0.26 ns _1 to within only 1%. We conclude that non- 
exponential decay of a single fluorescent nanosphere is a 
characteristic of the LDOS fluctuations in random pho- 
tonic media, and attribute the non-exponential decay to 
the str ong dependence of emission rate on dipole orien- 
tation [13| . 
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FIG. 2: (color) (a) Time-resolved fluorescence from a sin- 
gle fluorescent nanosphere embedded in PS random photonic 
medium with a thickness of 7.4(5) /im. The slow (purple) 
and the fast (blue) decay curves were measured on two dif- 
ferent single fluorescent nanospheres at the same depth d = 
3.3 /tm. The white lines through the data are lognormal fits. 
The most-frequent decay rate y m f and the width of the dis- 
tribution A7 are shown, (b) Time-resolved fluorescence of a 
single fluorescent nanosphere in a reference polymer sample 
(green) decays exponentially with a rate of 0.26 ns -1 . 

We quantitatively analyze the non-exponential time- 
resolved fluorescence curves with a distribution of emis- 
sion rates, using the method that successfully describes 
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time-resolved emission in ordered photonic media Q- 
By choosing a lognormal distribution of emission rates 
we are ensured of physical positive rates with only one 
extra adjustable parameter as compared to the single- 
exponential model. The parameters of the distribu- 
tion are the most-frequent emission rate J m f, which is 
the peak of the distribution, and A7, which is its 1/e 
width. The mean emission rate 7^ is equal to 7^ — 

7 m/ exp^3/4- [sinh -1 (A7/27 m/ )] 2 ^ 7]. From the anal- 
ysis of the emission curves of single probes inside PS sam- 
ples in Fig. [5J the slowest emission curve is described by 
7 M = 0.28 ns _1 while the fastest one has 7^ = 0.33 ns~L 
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FIG. 3: (color) Distribution of the emission rates for mate- 
rials with different scattering strengths (1/kl). Histograms 
of the emission rates from fluorescent nanospheres embedded 
in three different materials: (a) ZnO, (b) PS, and (c) non- 
scattering polymer layer. The width of the distribution 2a is 
indicated by the dotted lines. The full lines are Gaussian fits 
to the histograms. 

Figs. OHa) and (b) show histograms of emission rates 
for ZnO and PS random photonic media, respectively. 
As a measure of the width of the distributions we show 
twice the square root of the variance 2a — 2^/Var(7 Al ). 
The histogram for a non-scattering reference medium is 
shown in Fig. f5X c ) - For all our samples, we verified 
that both the mean and the variance of the emission rate 
are independent of sample thickness. The distribution of 
emission rates is very narrow in the reference samples, 
showing that the fluorescent spheres have negligible in- 
trinsic fluctuations. Interestingly, in moderately scatter- 
ing PS samples (1/kl =0.06), the width of distribution 
2cr = 0.05(1) ns _1 is already much larger than in the 
reference samples. In strongly scattering ZnO samples 
(1/kl =0.12), the distribution is twice as broad, 2a = 
0.10(2) ns~L We note in passing that the maximum of 
the distribution also shifts, probably due to a change in 
effective refractive index 2l|. The increase of the width 



tracks the increase of the scattering strength, in qualita- 
tive agreement with the theoretical prediction that LDOS 
fluctuations increase with the scattering strength [ToLfTl| . 

The main result of our work shown in Fig. 0] is the 
measurement of the relative variance of the emission rate 
ct 2 /(7 m ) 2 versus scattering strength (1/kl). For non- 



0) the relative variance van- 
For PS {1/ki = 0.06) and 



scattering samples (l/k£ - 
ishes (a 2 /(^) 2 < I0- 4 ). 
ZnO (1/kl — 0.12) samples, we find a relative variance of 
0.007(2) and 0.020(4), respectively. The relative variance 
clearly increases with the scattering strength. Strikingly, 
the LDOS fluctuations predicted from the point-like scat- 
terer model 
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are about a factor twenty too high com- 
pared to our measurements. We attribute this difference 
to the fact that scalar wave and point scatterer approx- 
imations were used in Ref. ll|, whereas in experiments 



the scatterers have finite size and the light waves have 
vector character. The vector character of light could have 
a strong effect as shown theoretically in Ref. 13j. 



Polymer PS ZnO 



o 

CD CO 
O CD 

1 .1 

'*= F 

CD 

en 



-^0.10 



0.05 



0.00* 




<> 



° CD 
CD 

CO CD 
-+— * 

CD 
> 



CD 
> 



c 
O 

CO 

Id E 

CD 



0.00 0.04 0.08 0.12 

Vki 



FIG. 4: (color) (a) Relative variance of the mean decay rate 
<7 2 /{7 M } 2 versus scattering strength 1/kl. The measured val- 
ues are shown as filled circles with error bars indicating the 
statistical uncertainty for N = 48 measurements. The empty 
diamonds are the upper bound to Co determined from the 
fluctuations of the emitted intensity. The full and dotted lines 
show the LDOS variance for point-like scatterers [llj and for 
the single-sphere-scatterer model, respectively. 



As the point-like scatterer model does not provide an 
accurate description of the immediate surroundings of the 
emitter, we propose an alternative model to estimate the 
fluctuations of the emission rate, in which only one scat- 
terer closest to the fluorescent nanosphere is considered. 
This is reasonable since in a dried sample, fluorescent 
nanospheres are always in contact with the surface of at 
least one scatterer. The calculation is done by adapt- 
ing the analytical model for the spontaneous emission of 
an electric dipole located near a spherical surface [22j ]. 
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We calculate the emission rate of a dye molecule at all 
positions inside a fluorescent nanosphere in contact with 
the surface of a PS or ZnO spherical scatterer for all 
possible orientations of the molecule. The relative vari- 
ance of the emission rate is obtained by avera gin g over 
the size distribution of the spherical scatterers [18|. The 
calculated values are 0.010(3) and 0.022(5) for PS and 
ZnO samples, respectively, in good agreement with the 
experimental data (Fig. SJ) (2^] - This result confirms the 
hypothesis that the nearest scatterer dominates LDOS 
fluctuations in random photonic media 1^, U | ■ 

It has been predicted that in infinite random media 
LDOS fluctuations are equivalent to Co fluctuations of 
emission intensity llj. To probe Co, we have measured 
the relative variance of the emission intensity of 20 indi- 
vidual fluorescent nanospheres inside 16.5 /^m-thick PS 
and 13.3 /urn-thick ZnO random photonic media. We 
note that the variance of the intensity provides an up- 
per bound to Co since it consists of both residual speckle 
variance and the desired Co contribution |18j ] . The inten- 
sity variance is therefore necessarily larger than Co- The 
relative variance of emission intensity is shown in Fig. 
[4j it is indeed always higher than the relative variance of 
the emission rates. This observation means that our data 
are consistent with the hypothesis that the LDOS fluctu- 
ations and Co fluctuations are equivalent in our samples. 

In conclusion, we observed fluctuations of the LDOS 
inside random photonic media through the spontaneous 
emission of individual fluorescent nanospheres embedded 
deep inside. Recently, an interesting study of emission 
rates near the interface of random materials has appeared 
25] . The relative variance of LDOS fluctuations increases 
with scattering strength and agrees with a theoretical 
model. We suggest that Co may also be responsible for 
the broad distribution of emission rates observed in pho- 
tonic crystals pj. In view of the overestimation of fluc- 
tuations by the point-like scatterer model, we anticipate 
that our observations will stimulate new studies of the 
LDOS in random and partly ordered photonic media. 
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